AD-A27S  837 


TECH.  MEMO 
AL.ivO  i.  I A  7 


UNLIMITED 


TECH.  M 
AERO  21 


’<>•!  VA-Av  A'ri 


ROYAL  AEROSPACE  ESTABLISHMENT 


A  STUDY  OF  FLOWS  OVER  HIGHLY-SWEPT  WINGS  DESIGNED 
FOR  MANOEUVRE  AT  SUPERSONIC  SPEEDS 

by 

P.  R.  Ashill 
J.  L.  Fulkcr 
M.  J.  Simmons 


October  1988 


rLECTF. 
JAN  19  B9D 


\%Sr- 


Procurement  Exocutivo,  Ministry  of  Defence 
Farnhorough,  Hampshire 


BEST 

AVAILABLE  COPY 


9  0  -01  18  0  21 

UNI.  I  MI  rF.P 


0056036 


CONDITIONS  OF  RELEASE 


BR-1 12012 


COPYRIGHT  (C) 
1988 

CONTROLLER 
HMSO  LONDON 


. . ***** .  Y 

Reports  quoted  are  not  necessarily  available  to  members  ol  the  public  or  to  commercial 
organisations. 


UNLIMITED 


ROYAL  AEROSPACE  ESTABLISHMENT 
Technical  Memorandum  Aero  2147 
Received  for  printing  21  October  1988 

A  STUDY  OF  FLOWS  OVER  HIGHLY-SWEPT  WINGS  DESIGNED 
FOR  MANOEUVRE  AT  SUPERSONIC  SPEEDS 

by 

P.  R.  Ashill 
J.  L.  Fulker 
M.  J.  Simmons 


SUMMARY 

A  wind-tunnel  investigation  into  supersonic  free-stream  flows  over  two 
wing-body  configurations,  having  wings  of  different  design,  suitable  for  combat 
aircraft,  is  described.  Both  wings  have  the  same  quasi-delta  planform  of  60 
inboard  leading-edge  sweep  and  the  same  4Z  thickness  distribution  but  have 
differing  camber  distributions.  Following  a  description  of  the  design  of  the 
wings,  the  test  procedures  are  discussed  and  the  general  features  of  the  flows 
at  conditions  close  to  those  for  sustained  manoeuvre  are  identified  and 
contrasted.  Comparisons  between  calculations  by  CFD  methods  and  measurement 
are  presented,  and  it  is  shown  that  a  multiblock  method  for  solving  the  Euler 
equations  is  suitable  for  designing  wings  for  efficient  manoeuvre  at  supersonic 
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INTRODUCTION 
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Future  combat  aircraft  will  be  required  to  manoeuvre  for  sustained  periods 
at  supersonic  speeds.  At  a  typical  manoeuvre  condition,  the  lift-dependent 
drag  is  a  major  contribution  to  aircraft  drag,  so  that  measures  taiten  to 
minimise  this  drag  component  may  confer  combat  superiority  on  the  aircraft. 

At  the  high  angles  of  incidence  needed  for  manoeuvre,  there  may  be 
regions  of  adverse  pressure  gradient  on  the  wing,  possibly  provoking  flow 
separation.  Such  separations  can  either  be  exploited,  as  with  'vortex 
flaps',  or  be  eliminated  by  careful  design.  In  either  case,  the  aim  is  to 
maximise  the  thrust  from  suctions  on  forward-facing  surfaces  of  the  wing 
and  thus  to  minimise  drag. 

This  paper  Is  concerned  with  wings  designed  for  attached  flow, 
although  the  implications  for  drag  of  Halted  regions  of  separation  are 
considered.  Some  principles  underlying  the  design  of  conical  wings  with 
attached  flow  have  been  described  by  Miller  et  cl  (Ref  1)  and  are 
Illustrated  In  Fig  1  by  sketches  of  a  wing  spanwlse  pressure  distribution 
and  conical  or  'cross-flow'  streamlines,  ie  the  intersections  of  the  stream 
surfaces  with  a  sphere  centred  on  the  wing  apex  (Ref  2).  In  the  type  of 
flow  described  by  Miller  et  at,  the  component  of  free-stream  flow  normal  to 
the  leading  edge  Is  subsonic  but,  at  high  Incidence,  the  cross  flow  accel¬ 
erates  around  the  leading  edge  to  become  supercritical  on  the  upper  sur¬ 
face.  The  wing  is  cambered  In  the  spanwlse  sense  to  reduce  both  the  shock 
strength  dose  to  the  wing  surface  (Ref  1)  end  the  penetration  of  the  shock 
lnco  the  flovfield,  thereby  achieving  the  dual  aim  of  minimising  or  elimin¬ 
ating  the  wave  drag  due  to  cross-flow  shocks  and  suppressing  the  associ¬ 
ated  separations.  Camber  also  Increases  the  forward-feeing  area  available 
to  provide  thruat.  Unsustainable  suction  peaks  at  the  leading  edge,  both 
on  and  off  design,  are  prevented  by  the  use  of  rounded  leading  edges. 
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Although  these  ideas  are  applicable  to  conical  vlngs,  they  are  useful 
for  interpreting  flows  over  lifting  wings  of  quasi-delta  planform  and  with 
sharp  trailing  edges,  since,  as  will  be  seen  later,  the  flows  on  the  upper 
surface  of  these  wings  are  of  a  conical  nature. 

Since  the  velocities  on  the  upper  surface  can  be  large  at  manoeuvre 
conditions,  methods  based  on  linearised  theory  are  unsuitable  for  the 
design  of  efficient  manoeuvre  wings.  Methods  which  solve  the  exact  poten¬ 
tial  equation  are  valid  for  flows  with  weak  shocks  and,  as  such,  have  been 
used  to  design  both  conical  and  non-conical  wings  for  shock-free  flows  on 
the  upper  surface  (Refs  1  and  3).  However,  these  methods  do  not  represent 
boundary-layer  displacement  effects  and  neglect  the  increase  in  entropy 
through  the  shock:  thus  potential-flow  methods  may  be  seriously  In  error 
when  cross-flow  shocks  of  significant  strength  are  present  and  hence  may 
limit  the  freedom  to  make  design  compromises  or  to  study  the  effects  of 
departures  from  the  design  condition.  Methods  are  available  for  solving 
the  Euler  equations  (Refs  4  and  5)  which  represent  the  changes  In  flow 
quantities  across  shocks,  although,  as  with  the  potential-flow  methods, 
they  do  not  represent  boundary-layer  effects. 

Therefore  detailed  experimental  studies  are  needed  of  flows  over 
supersonic  manoeuvre  wings  to  validate  Computational  Fluid  Dynamic  (CFD) 
methods  and  to  identify  flow  features  not  adequately  modelled.  This  paper 
describes  an  Investigation  with  these  alms,  which  was  performed  In  the 
8ft  x  8ft  Wind  Tunnel,  RAE  Bedford,  on  half  models  of  two  wing-body  con¬ 
figurations,  each  having  the  same  quasi-delta  planform  of  60°  inboard 
leading-edge  sweep  and  thickness  form  of  41  thickness/chord  ratio  but  with 
differing  camber.*  The  large  site  of  the  models  enabled  them  to  be  manu¬ 
factured  to  the  desired  accuracy  and  also  allowed  flow  measurements  to  be 
made  In  considerable  detail.  Manufacturing  accuracy  Is  particularly 
Important  in  the  highly-curved  region  of  the  leading  edge  which  controls 
the  development  of  the  flow  on  the  upper  surface. 

After  briefly  describing  the  methods  used  to  design  the  two  wings  In 
section  2,  the  paper  continues  with  a  description  of  the  wind-tunnel  tests 
In  section  3.  A  general  description  of  the  flows  Is  provided  In  section  4, 
and  comparisons  between  predictions  by  CFD  methods  and  measurement  are 
presented  In  section  S.  Section  S  concludes  with  an  assessment  of  the  drag 
characteristics  of  the  two  wings  and  considers  the  Implications  for  the 
aerodynamic  design  of  wings  for  supersonic  combat  aircraft. 


2  WXgG  DESIGN 

The  two  wings  are  Illustrated  In  Fig  2;  Wing  A  has  a  complex  camber  surface 
with  camber  in  both  spanwlse  and  atreamwlse  directions  while  Wing  B  is  of 
conical  camber  with  the  apex  at  the  leading  edge  of  the  wing-body  junction. 
The  thickness  distribution  of  stresmwlse  sections  is  a  two-parameter  type 
with  maximum  thickness /chord  ratio  4Z,  as  noted  before.  The  other 


*  Recently,  tests  were  completed  on  a  third  member  of  the  family  of  wings 
having  a  symmetrical  section.  Analysis  of  these  tests  Is  not  complete 
and  so  no  further  reference  Is  made  to  this  wing  except  In  section  5 
where  data  from  It  are  used  in  the  analysis  of  the  lift-dependent  drag  of 
the  other  two  wings. 
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parameter  was  chosen  to  vary  across  the  wing  with  the  aim  of  keeping  the 
maximum-thickness  locus  well  swept  and  having  a  relatlvely-large  nose 
radius  on  the  outer  wing  to  avoid  large  suction  peaks  there. 

Fig  3  shows  the  spanwlse  variation  of  the  ratio  of  the  section  nose 
radius  pl  to  local  chord  c  .  Inboard  of  mid  span,  the  nose-radius  ratio 
Is  of  similar  magnitude  to  those  of  transonic  aerofoils  with  the  same 
thickness.  Further  outboard,  the  nose-radius  ratio  increases  rapidly  with 
spanwlse  distance,  varying  between  50Z  and  80Z  span  Is  a  way  similar  to  the 
trend  for  a  conical  wing  for  which  the  leading-edge  suction  peak  is 
constant  across  the  span  (Fig  3). 

Also  Illustrated  In  Fig  3  Is  the  distribution  of  nose  radius  of  a 
delta  wing  designed  by  Hood  and  Bauer  (Ref  6).  This  design  has  a  maximum¬ 
thickness  locus  which  Is  parallel  to  the  leading  edge  (cf  the  maximum- 
thickness  locus  of  the  present  wing.  Fig  3)  and  uses  sections  with  non-zero 
base  thickness  outboard  of  66Z  span  to  derive  a  larger  nose  radius  In  this 
region.  This  spanwlse  variation  of  nose  radius  Is  similar  to  that  of  the 
present  wing.  Using  a  calculation  method  for  solving  the  full-potential 
equation.  Hood  and  Bauer  found  that  their  wing  has  lower  lift-dependent 
drag  than  the  corresponding  wing  with  a  constant  section  across  the  span. 
They  argued  that  this  comes  about  because  the  former  wing  has  a  larger 
forward-facing  area  than  the  latter  for  approximately  the  same  conical-type 
pressure  distribution. 

Wing  A  has  a  camber  surface  based  on  that  of  an  earlier  design  which 
was  studied  as  part  of  an  extensive  programme  of  research  conducted  at  RAE 
In  the  1950's  on  slender  wings  with  sharp  leading  edges  suitable  for  super¬ 
sonic  transport  aircraft.  Designed  using  linearised  theory,  (Ref  7)  the 
original  wing  was  of  71*  leading-edge  sweep  and  was  cambered  so  that,  at 
the  cruise  Mach  number  2.2,  the  leading  edge  Is  an  attachment  line  at  a 
lift  coefficient  Cl  "  0.05.  Chord  and  span  loadings  were  chosen  with  the 
dual  aim  of  achieving  low  lift-dependent  drag  and  of  trimming  the  aircraft 
at  the  cruise  lift  coefficient  0.1.  Wing  A  retains  the  spanwlse  camber 
distribution  of  the  original  wing  but  with  the  ordinates  scaled  to  give  a 
leading-edge  sweep  of  60°.  According  to  the  Prandtl-Glauert  similarity  law 
this  Implies  that  the  leading  edge  of  Wing  A  Is  an  attachment  line  at 
Cl  3  0.1  and  at  a  Mach  number  of  1.5.  This  Mach  number  Is  typical  of 
supersonic  sustained-turn  requirements  but  the  lift  is  about  one  third  of 
that  needed  at  the  tropopause.  However,  as  Is  shown  later,  the  rounded 
leading  edge  ensures  that  the  flow  remains  attsched  there  for  lift  coef¬ 
ficients  up  to  and  beyond  those  of  Interest. 

The  camber  surface  of  Wing  B  was  designed  using  the  COREL  code  which 
solves  numerically  the  exact  potential  equation  for  conical  flow  (Ref  8) 
and  Includes  an  approximate  allowance  for  non-conlcal  effects  based  on  the 
theory  of  plane  waves  (Ref  9).  Developed  at  Grumman  Aerospace  and  made 
available  to  RAE  by  NASA  under  the  terms  of  an  MOD/NASA  collaboration,  the 
method  Is  straightforward  to  use,  with  a  number  of  useful  features  for 
design  purposes.  In  Its  present  form,  the  method  cannot  represent  a  non- 
conlcal  body,  and  consequently  the  'net-wing'  spproxlmatlon  was  used  (le 
the  wing  plane  of  symmetry  was  taken  at  the  body  side). 

Fig  4  illustrates  spanwlse  pressure  distributions  calculated  by  this 
method  for  both  wings  at  the  axial  station  x/co  «  0.705  (as  defined  In 
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Pig  5)  and  for  M  •  1.6.  Here  the  static-pressure  coefficient  Cp  is 
plotted  against  spanvlse  distance  from  the  body  centre-line  divided  by 
local  semi-span,  ETA.  On  the  left-hand  aide  of  the  figure  are  shown 
pressure  distributions  corresponding  to  the  design  condition  of  Wing  B,  le 
at  a  lift  coefficient  based  on  exposed  wing  area  -  0.4. 

At  x/c0  »  0.705,  the  flow  over  Wing  B  at  the  'design'  condition  is 
predicted  to  be  shock-free  on  the  upper  surface  while  a  strong  shock  Is 
Indicated  on  the  upper  surface  of  Wing  A.  For  lower  lift  coefficients  or 
further  upstream.  Wing  B  is  predicted  to  have  a  re-expansion  on  the  upper 
surface  Inboard  of  the  Initial  compression,  followed  further  Inboard  by  a 
more  severe  compression.  However,  calculations  by  a  three-dimensional 
boundary-layer  method  (Ref  10)  suggest  that  separation  would  be  avoided  on 
the  upper  surface  for  stations  as  far  aft  as  x/c0  -  0.705  and  for  lift 
coefficients  up  to  0.4  at  M  •  1.6. 

Fig  4  also  shows  calculated  pressure  distributions  for  C^g  ”0.1 
corresponding  to  that  for  supersonic  dash.  In  the  design  of  Wing  B, 
special  efforts  were  made  to  prevent  separation  on  the  lower  surface  at 
this  condition,  and,  in  this  respect,  Fig  4  shows  that  the  adverse  pressure 
gradients  on  this  surface  are  smaller  on  Wing  B  than  on  Wing  A.  The  extent 
to  which  the  design  alms  of  Wing  B  have  been  met  Is  discussed  in  section  4. 


3  MODELS  AND  TEST  TECHNIQUES 

The  layout  of  both  models  in  the  wind  tunnel  Is  Illustrated  In  Fig  5;  each 
wing  was  mounted  low  on  a  half  body  which  Is  cylindrical  and  nearly  rec¬ 
tangular  In  cross  section  where  the  wing  Intersects  the  body.  Further 
upstream,  the  body  becomes  of  circular  cross  section  as  the  nose  Is 
approached. 

In  each  case,  the  model  was  attached  to  a  strain-gauge  balance  outside 
the  working  section  for  the  measurement  of  overall  forces.  The  centre-line 
of  the  body  was  offset  approximately  25  in  from  the  tunnel  sidewall  to 
minimise  Interference  between  the  model  and  the  sidewall  boundary  layer. 
Inspection  of  the  wings  revealed  a  high  standard  of  surface  finish 
(0. 1-0.2  urn)  and  that  manufacturing  errors  In  the  ordinates  are  small  being 
typically  of  the  order  0,1  in  (0.00006  Cq). 

As  shown  In  Fig  5,  pressure  tappings  were  provided  at  nine  streaasrlse 
stations  on  each  wing  and  on  the  side  of  the  body  along  Its  centre-line. 
Model  pressures  were  measured  using  nine  transducers  each  Installed  In  a 
48-way  pressure  switch  and  having  a  working  range  of  ±1.7  bar.  Transducer 
calibration  was  effected  during  each  scan,  using  tunnel  total  and  static 
pressures  as  known  datum  conditions. 

The  tunnel  has  good  control  of  stagnation  pressure,  temperature  and 
humidity  down  to  a  frost  point  of  233  K.  The  various  tunnel  and  reference 
pressures  were  measured  by  Individual  self-balancing  capsule  manometers 
having  highly  stable  calibrations  and  zeros.  The  manometers  recording 
stagnation  pressure  have  a  resolution  of  0.34  mbar  and  those  monitoring 
wall  static  pressures  0.17  mbar.  Tunnel  total  temperature  was  controlled 
such  that  excursions  In  balance  temperature  during  a  test  were  minimised. 
This  precaution  together  with  corrections  To  the  balance  zeros  for  change 
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in  temperature  ensured  a  high  standard  of  repeatability  of  balance  force 
measurements.  For  force  and  pressure  measurement,  frost  point  was  main¬ 
tained  below  243  K. 

The  correction  to  nominal  Mach  number  was  obtained  directly  from  che 
tunnel  calibration  and  for  the  conditions  studied  Increases  Mach  number  by 
0.005. 

In  Ref  11  it  is  shown  that  measurements  of  both  forces  and  pressures 
were  made  with  a  high  standard  of  accuracy  and  repeatability.  Measurements 
of  flow  pitch  angle  indicate  that  the  overall  error  in  model  incidence  is 
small  and  probably  in  the  range  ±0.02°. 

Boundary-layer  transition  was  fixed  by  narrow  bands  of  sparsely- 
distributed  ballotlnl,  cemented  to  the  model  by  epoxy  resin.  On  both  sur¬ 
faces  of  the  wing  the  band  was  25.4  mm  from  the  leading  edge  in  plan  view 
and  was  5.08  mm  wide.  On  the  body  the  band  was  12.7  mm  wide  and  was 
located  101.6  mm  along  the  surface  downstream  of  the  nose.  In  all  cases, 
the  ballotlnl  were  within  the  range  of  diameters  0.21-0.25  mm. 

Flow  visualisations  were  performed  using  surface  oil  flows  and  vapour 
screens.  The  oil  was  supplied  to  the  upper  surface  of  each  wing  through 
holes  drilled  normal  to  the  surface  in  three  tubes  laid  across  the  span 
(Fig  5).  Photographs  of  the  oil  flows  were  taken  through  a  window  in  the 
working-section  roof  to  the  rear  of  the  model.  The  vapour  screens  were 
illuminated  by  a  vertical  fan  of  laser  light  produced  by  reflecting  a  beam 
from  a  2  H  continuous-wave,  Argon  laser  with  an  oscillating  mirror.  The 
streanwlse  position  of  the  fan  could  be  changed  by  rotating  the  mirror 
about  a  vertical  axis;  thus  the  fan  was  normal  to  the  tunnel  axis  at  only 
one  streamwlse  station  (see  later).  The  development  of  the  system  and  the 
procedures  used  to  obtain  satisfactory  vapour  screens  are  described  in 
Ref  12.  Photographs  of  the  screens  were  taken  by  cameras  at  three  pos¬ 
itions;  (a)  in  a  pod  approximately  1.2  Cq  downstream  of  the  wing  trailing 
edge,  (b)  in  che  Schlleren  cavity  opposite  the  model  and  (c)  in  the 
working-section  roof  to  the  rear  of  che  model.  A  typical  photograph  by  the 
camera  in  che  Schlleren  cavity  is  shown  in  Fig  6.  Using  arguments  put  for¬ 
ward  convincingly  by  McGregor  (Ref  13),  it  is  possible  to  identify  the  main 
features  of  che  flow  from  such  a  vapour  screen  as  indicated  in  Fig  6. 

The  tests  were  performed  at  a  unit  Reynolds  number  of  13.1  x  10®/m  (le 
at  a  Reynolds  number  based  on  geometric  chord  of  12.7  x  10®)  and  for  the 
Mach  numbers  1.4,  1.6  and  1.8. 

Tests  were  also  made  on  the  body  alone  for  the  same  range  of  con¬ 
ditions  as  the  wing-body  combinations,  with  the  aim  of  providing  a  basis 
for  analysing  the  overall  forces  on  the  wings. 

4  GENERAL  DESCRIPTION  OF  FLOWS 

This  section  describes  and  contrasts  the  flows  over  the  two  wings  using 
measurements  of  wing  pressures  and  flow  visualisations.  The  discussion  is 
confined  to  M  «  1.6,  since  the  flows  observed  at  this  Mach  number  are  simi¬ 
lar  In  character  to  those  found  at  the  other  Mach  numbers  studied. 
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Fig  7  Illustrates  wing  pressure  distributions  at  M  ■  1.6  for  a  number  of 
Incidences,  a  ,  and  at  x/c0  -  0.705.  Close  to  zero  lift  there  Is  a 
pronounced  suction  peak  near  the  leading  edge  on  the  lower  surface.  As 
Incidence  Increases,  this  peak  diminishes  -  presumably  because  the  attach¬ 
ment  line  moves  from  the  upper  to  the  lower  surface  -  and  the  suctions  on 
the  upper  surface  Increase.  The  upper-surface  suctions  Increase  gradually 
with  distance  from  the  leading  edge  until  pressure  recovery  begins  at 
between  70  and  80Z  span,  the  pressure  recovery  being  effected  partly  by  a 
shock  for  Incidences  above  about  5*.  At  all  incidences,  the  pressure 
distributions  suggest  that  the  flow  Is  attached  outboard  of  the  shock  on 
the  upper  surface.  However,  at  lift  coefficients  above  about  0.2,  separ¬ 
ation  occurs  at  the  foot  of  the  shock  as  Indicated  by  the  characteristic 
bulge  in  the  pressure  distribution. 

Figs  8,  9  and  10  Illustrate  Isobar  patterns  on  the  upper  surface 
(plotted  at  Intervals  of  0.05  In  Cp  unless  otherwise  Indicated)  along 
with  sketches  from  photographs  of  oll-fllament  flows  and  vapour  screens  for 
o  -  5.3*  (CL  -  0.2),  7.8*  (CL  -  0.3)  and  10.3*  (CL  -  0.4).  As  noted  In  the 
Introduction,  these  flows  are  of  conical  character,  eg  the  upper-surface 
shock  is  closely  aligned  with  a  generator  of  the  'gross'  wing  except  within 
the  aft  20Z  chord  at  the  two  highest  Incidences. 

In  the  vapour  screens  of  Figs  8  to  10  and  In  subsequent  figures,  the 
bow  wave  is  indicated  by  a  line  outboard  of  the  wing  and  roughly  normal  to 
the  edges  of  the  laser  fan. 

At  a  •  5.3“  the  Isobars  outboard  of  the  shock  are  generally  well  swept 
and  there  Is  no  evidence  In  the  flow  visualisations  of  leading-edge  separ¬ 
ation.  A  shock  above  the  upper  surface  Is  Indicated  by  the  pronounced 
turning  of  the  oil  filaments,  the  concentration  of  Isobars  and  the  lines 
roughly  normal  to  the  wing  surface  in  the  vapour  screens  at  approximately 
70Z  span.  Shock-Induced  separation  Is  less  certain,  although  a  thin  dark 
region  In  the  vapour  screen  just  above  the  wing  suggests  that  this  Is  a 
possibility. 

At  u  ■  7.8*,  (CL  »  0.3)  the  Isobars  near  the  leading  edge  converge 
with  increasing  axial  distance.  However,  the  vapour  screens  do  not  point 
to  leading-edge  separation  anywhere  on  the  wing.  The  deflection  of  the  oil 
filaments  through  the  shock  Is  greater  than  at  a  •  5.3*,  and  a  definite 
Indication  of  shock-induced  separation  Is  provided  by  the  outboard  deflec¬ 
tion  of  the  oil  filaments  Inboard  of  the  shock.  This  Inference  Is  con¬ 
firmed  by  the  vapour  screens  which  each  reveal  a  dark  bulbous  region 
Inboard  of  the  shock  Indicative  of  a  separation  bubble.  For  stations  aft 
of  about  60Z  chord,  the  termination  of  this  separation  or  'reattachment'  Is 
Indicated  by  the  clustering  together  of  Isobars  at  about  502  span.  A 
smaller  dark  region  In  the  vapour  screens  Inboard  of  the  main  separation 
appears  to  have  Its  origin  at  the  junction  of  the  wing  leading-edge  and  the 
body,  suggesting  that  this  Is  a  separation  resulting  from  wing-body  Inter¬ 
ference.  However,  this  feature  does  not  have  a  significant  effect  on  wing 
pressures.  The  vapour  screens  also  show  that,  at  a  given  axial  station  and 
Just  above  the  wing,  the  Inclination  of  the  shock  with  respect  to  the  wing 
surface  is  smaller  at  i  ■  7.8*  than  at  5.3*.  Furthermore,  the  deflec¬ 
tion  of  the  flow  normal  to  the  wing  surface  as  it  passes  through  the  shock 
Increases  noticeably  with  Increase  In  Incidence.  The  theory  of  oblique 
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shocks  (Ref  14)  shows  that  these  observations  are  consistent  with  the  shock 
tt  «■  7.8°  being  stronger  than  that  for  a  •  5.3®  as  would  be  expected. 

The  loss  in  isobar  sweep  near  the  leading  edge  aft  of  803  chord  is 
particularly  pronounced  at  a  »  10.3°  (Cl  “  0.4)  but  again  there  is  no  sign 
of  leading-edge  separation.  A  further  reduction  in  shock  inclination  just 
above  the  wing  is  seen  in  the  vapour  screens  along  with  an  Increase  in  both 
the  deflection  of  the  flow  normal  to  the  wing  surface  after  passing  through 
the  shock  and  the  cross-sectional  area  of  the  separation.  In  addition,  a 
second  limb  is  apparent  in  the  shock  above  the  main  separation  region  in 
the  furthest-aft  vapour  screen,  and,  in  all  the  vapour  screens,  dark 
patches  are  evident  just  above  the  inboard  extremity  of  the  main  separation 
indicating  regions  of  re-expanslon  followed  by  weak  shocks.  The  concen¬ 
tration  of  Isobars  Inboard  of  the  shock,  symptomatic  of  reattachment, 
extends  further  forward  than  for  a  ■  7.8®. 

Near  the  intersection  of  the  upper-surface  shock  with  the  trailing 
edge,  the  oil  filaments  show  that  the  surface  flow  has  a  pronounced  span- 
wise  component.  This  effect  is  probably  due  to  the  combined  influence  of 
the  tralllng-edge  and  upper  surface  shocks  on  the  boundary  layer. 

The  features  noted  above  become  more  obvious  with  further  Increase  of 
incidence  and  at  values  of  incidence  greater  than  about  12°  the  main 
separation  clearly  becomes  detached  from  the  wing  surface  and  appears  to  be 
of  the  vortex-sheet  type. 

4.2  Wing  B 

Fig  11,  which  corresponds  to  Fig  5  for  Wing  A,  illustrates  the  main 
features  of  the  pressure  distributions  of  Wing  B  at  x/c„  ■  0.705.  As  with 
Wing  A,  a  suction  peak  is  evident  near  the  leading  edge  on  the  lower  sur¬ 
face  at  low  lift.  However,  as  Intended,  this  peak  is  of  lower  magnitude 
than  that  of  the  corresponding  peak  on  Wing  A.  This  aspect  is  considered 
again  in  section  5. 

The  upper-surface  pressure  distributions  of  Wing  B  also  differ  in 
character  from  those  of  Wing  A.  By  design.  Wing  8  has  higher  suctions  on 
the  outboard  5-103  span  than  does  Wing  A,  with  some  lsentroplc 
recompression  evident  at  the  three  largest  incidences.  Since  these  suc¬ 
tions  act  on  a  forward-facing  surface,  they  are  a  desirable  feature  of  the 
design  of  Wing  B.  However,  an  unintended  re-expanslon  occurs  further 
Inboard  leading  to  a  shock  forming  for  a  ■  5.4®  and  higher  Incidences. 
Reasons  for  this  are  discussed  in  section  5. 

Figs  12,  13  and  14  show  Isobars  and  flow  visualisations  for  a  ■  5.4®, 
7.5®  and  10®  (Cl  *  0.2,  0.3  and  0.4).  As  with  Wing  A,  the  upper-surface 
flows  have  conical  features  but,  in  this  ease,  the  upper-surface  shock  lies 
close  to  a  generator  of  the  'net'  wing.  The  flow  visualisations  of  the  two 
wings  are  broadly  similar;  however.  Wing  B  does  not  exhibit  the  dark  patch 
Inboard  of  the  main  separation  noticed  in  the  vapour  screens  of  Wing  A  and, 
furthermore,  Wing  B  is  less  affected  by  separation  than  Wing  A,  as  may  be 
seen  by  comparing  Fig  9  with  Fig  13  and  Fig  10  with  Fig  14.  These  compari¬ 
sons  reveal  that  (a)  the  shock-induced  separation  bubble  of  Wing  A  is 
bluffer  than  that  of  Wing  B  and  (b)  the  inclination  of  the  upper-surface 
shock  to  the  wing  surface  is  smaller  for  Wing  A  chan  for  Wing  B,  in  both 
cases  for  a  given  axial  position  and  lift  coefficient.  The  lower  shock 


r 


Inclination  of  Wing  A  compared  to  that  of  Wing  B  Is  consistent  with  the 
former  wing  having  a  stronger  upper-surface  shock  than  the  latter  and  this 
Is  confirmed  by  calculations  of  shock  strength  from  the  measured  pressure 
distributions  to  be  described  belov. 

Comparison  between  Figs  9  and  14  shows  that  Wing  B  has  roughly  the 
same  size  and  shape  of  upper-surface  separation  at  -  0.4  as  Wing  A  has 
at  Cl  ■  0.3.  A  more  precise  Indication  of  the  relative  significance  of  the 
separations  of  the  two  wings  may  be  obtained  by  studying  separation  'bluff¬ 
ness',  as  shown  by  the  vapour  screens.  Bluffness,  defined  as  the  ratio  of 
maximum  bubble  height  relative  to  the  wing  surface,  h  ,  to  bubble  width 
1  ,  Is  a  useful  measure  of  the  effect  that  the  separation  has  on  wing 
pressures.  Fig  15  shows  values  of  separation  bluffneas  Inferred  from  the 
vapour  screen  labelled  x/c0  -  0.789  (Fig  14)  for  a  range  of  angles  of  inci¬ 
dence.  The  plots  of  this  ratio  against  C^  show  that  Wing  B  has  a  lift 
coefficient  of  about  0.08  higher  than  that  of  Wing  A  for  a  given  value  of 
h/1  ,  confirming  the  remark  above.  Plotted  In  the  same  figure  against  M(j 
(the  Mach  number  of  the  flow  Just  upstream  of  and  normal  to  the  locus  of 
the  shock  on  the  wing  surface.  Inferred  from  the  measured  pressure  distri¬ 
butions  by  simple  sweep  theory  on  the  assumption  that  the  entropy  rise 
through  the  bow  shock  can  be  Ignored)  the  values  of  the  ratio  h/1  for  the 
two  wings  appear  to  collapse  onto  a  single  curve.  Intersecting  the  abscissa 
at  %  J  1.3.  In  other  words,  shock-induced  separation  is  only  significant 
for  values  of  %  »  1.3. 

In  one  small  part  of  the  flow.  Wing  B  shows  a  stronger  vlscous- 
lnvlscld  Interaction  than  does  Wing  A:  this  area  is  dose  to  the  Intersec¬ 
tion  of  the  upper-surface  shock  with  the  trailing  edge.  The  oil  flows  of 
Figs  9  and  13  reveal  that  the  spanwlse  flow  In  this  area  is  greater  for 
Wing  B  chan  for  Wing  A.  It  is  suggested  that  this  difference  Is  due  to  the 
tralllng-edge  shock  of  Wing  B  being  stronger  than  that  of  Wing  A,  an 
assertion  supported  by  the  measurements  of  pressures  near  the  trailing 
edges  of  the  two  wings. 


5  COMPARISONS  BETWEEN  CFD  METHODS  AND  MEASUREMENT 

In  this  section,  comparisons  are  presented  between  predictions  by  CFD 
methods  and  measurement  and  the  drag  characteristics  of  the  two  wings  are 
assessed.  The  methods  are  COREL,  already  referred  to  in  section  3,  and  the 
Euler-Multiblock  code,  referred  to  below  as  EM  for  brevity.  EM  was  devel¬ 
oped  jointly  by  Aircraft  Research  Association,  Bedford  (Ref  5)  and  British 
Aerospace,  Fllton  (Ref  15)  and  Is  based  on  the  finite-volume  method  for 
solving  the  Euler  equations  of  Jameson  et  al  (Ref  4). 

Boundary-layer  effects  are  not  represented  In  either  COREL  or  EM. 
However,  since  the  flows  considered  are  supersonic,  boundary-layer  effects 
are  expected  to  be  Important  only  In  the  region  of  the  shock. 

As  noted  in  section  3,  there  Is  no  provision  In  COREL  to  represent  a 
non-conlcal  body  and,  because  of  this.  Wing  3  was  designed  using  the  'net- 
wing'  approximation.  Figs  16  and  17  show  comparisons  of  calculations  by 
this  method  and  measurement  for  Wing  B  at  M  ■  1.6  and  a  ■  5.4*  and  7.5* 
(corresponding  to  measured  overall-lift  coefficients  of  0.2  and  0.3). 

These  figures  suggest  that  the  'net-wing'  approximation  Is  misleading  In 
the  way  that  It  falls  to  represent  adequately  the  re-expanalon  of  the  flow 
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on  Che  upper  surface.  Some  Improvement  in  agreement  between  COREL  and 
measurement  Is  obtained  If  Che  'gross-wing'  approximation  Is  used  (le  the 
wing  is  extended  within  the  body).  Therefore,  in  any  future  use  of  wing- 
alone  methods  to  design  highly-swept  wings  for  supersonic  manoeuvre,  it  Is 
recommended  that  the  'gross-wing'  approximation  be  used. 

A  further  development  of  COREL  In  which  the  non-conical  potential 
equation  Is  solved,  NCOREL,  (Ref  16)  enables  the  flow  over  wing-body  com¬ 
binations  to  be  calculated.  Computations  by  this  code  of  flows  over  Wing  A 
(with  the  body  represented)  have  been  made  by  NASA,  and  comparisons  between 
prediction  and  measurement  are  presented  in  Ref  17.  These  predictions  are 
comparable  In  accuracy  to  those  by  EM  given  here. 

The  topology  used  in  the  calculations  by  EM  described  In  this  paper 
has  'C'  and  'H'  structures  wrapped  around  streaowlse  and  spanwlse  sections 
of  the  configuration,  respectively,  and  an  'O'  structure  around  the  body 
cross-section.  The  number  of  blocks  used  is  42,  and  there  are  702  cells  on 
each  exposed  surface  of  the  wing.  All  the  solutions  are  numerically  well- 
converged,  and  in  Ref  11  it  is  shown  that  results  for  wing  pressure  distri¬ 
butions  are  not  sensitive  to  either  grid  topology  or  density  except 
possibly  close  to  the  wing  leading  edge. 

EM  predictions  and  measurement  of  spanwlse  pressure  distributions  for 
Wing  A  are  shown  in  Fig  18a-d  for  M  »  1.6  and  a  ■  0.29°,  3.31°,  5.32° 
and  7.80°  corresponding  to  measured  values  of  Cl  0,  0.1,  0.2  and  0.3.  At 
low  lift.  Cl  “  0,  (Fig  18a)  EM  underestimates  the  suction  peak  on  the  lower 
surface  except  at  x/c„  "  0.874  where  the  suction  peak  In  the  measurement 
appears  to  have  collapsed,  perhaps  owing  to  flow  separation.  At  higher 
lift.  Cl  “  0.1  and  0.2,  EM  gives  reasonable  predictions  of  the  pressures 
except  ac  Cl  “  0.2  Inboard  of  the  shock  on  the  upper  surface  where  the 
calculation  method  overestimates  the  suctions.  This  tendency  Is  more 
obvious  at  Cl  "  0.3  where  shock-induced  separation  Is  evident.  Undulations 
in  the  calculated  pressure  distributions  Inboard  of  the  shock  are  attri¬ 
buted  to  discretisation  errors  In  EM. 

Corresponding  comparisons  are  shown  In  Fig  19a-d  for  Wing  B  and  Indi¬ 
cate  thac  EM  predicts  the  pressure  distributions  of  this  wing  reasonably 
well  except  close  to  the  upper-surface  shock.  The  pressure  distributions 
of  Wing  A  differ  significantly  from  those  of  Wing  B,  as  noted  above,  and. 

In  both  cases,  EM  predicts  the  pressure  distributions  with  reasonable 
accuracy.  This  observation,  combined  with  a  knowledge  of  critical  con¬ 
ditions  for  shock-induced  separation  dlacussed  in  section  4,  points  to  the 
use  of  EM  as  a  design  tool.  A  successful  aspect  of  the  design  of  Wing  B  In 
avoiding  separation  on  the  lower  surface  at  low  lift  Is  Illustrated  by  the 
relatively  good  agreement  between  calculation  and  measurement  for  the 
pressures  o-.  this  surface  In  Fig  19a. 

Comparisons  between  predicted  and  measured  overall  forces  and  pitching 
moment  are  described  In  Ref  11.  In  this  analysis,  forces  and  pitching 
moment  on  the  body  alone  are  subtracted  from  those  of  the  wing-body  combi¬ 
nation  at  each  angle  of  Incidence,  thus  removing  extraneous  effects  such  as 
those  due  to  the  boundary-layer  diverter  between  the  body  and  the  tunnel 
sidewall.  The  results  of  this  procedure  are  called  notional  wing-alone 
forces.  An  analoguous  procedure  Is  used  Ip  the  calculations  by  EM,  along 
with  the  approximation  (of  slender-body  theory)  that  the  forces  on  the 
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cylindrical  pare  of  Che  body  alone  may  be  Ignored.  Generally,  the  agree¬ 
ment  between  calculation  and  measurement  of  lift  Is  reasonable  except  for 
flows  with  shock-induced  separation  for  which  EM  overestimates  the  lift  at 
a  given  Incidence.  As  already  shown,  this  discrepancy  arises  because,  for 
such  flows,  the  predicted  level  of  suctions  on  the  upper  surface  Inboard  of 
the  shock  are  higher  than  those  measured.  At  low  lift  (C^  <  0.1),  EM 
underestimates  drag  coefficient  by  about  0.001,  possibly  because  of  the 
failure  of  the  method  to  resolve  sufficiently  accurately  the  suction  peaks 
on  the  lower  surface*  or,  in  the  case  of  Wing  A,  because  of  lower-surface 
separation.  Pitching  moment  Is  not  predicted  accurately,  perhaps  because 
of  Che  known  sensitivity  of  pitching  moment  to  errors  in  wing  pressures  and 
also  because  of  the  neglect  in  the  calculation  method  of  forces  on  the 
cylindrical  part  of  the  isolated  body.  However,  the  difference  in  zero- 
lift  pitching  moment  coefficient  (based  on  geometric  mean  chord)  between 
Che  two  wings  (A-8)  is  estimated  with  reasonable  accuracy,  the  value  from 
experiment  being  »  0.0165  and  by  calculation  0.0145.  This  level  of 

accuracy  should  be  adequate  for  the  assessment  of  Che  relative  trim  drags 
of  different  wings. 

Fig  20,  taken  from  Kef  11,  shows  differences  in  drag  and  axial  force 
between  the  two  wings,  the  suffix  w  referring  to  notional  wing-alone 
forces.  At  low  lift  (C^  a  CLw  <  0.1),  the  predictions  of  drag  difference 
are  in  reasonable  accord  with  measurement  (Fig  21a).  However,  at  values  of 
CLw  above  0.2,  the  measured  difference  increases  relative  to  the  predicted 
values.  Fig  21b  reveals  that  the  trend  of  ACAw  with  CLw  is  similar  to 
that  measured,  implying  chat  the  discrepancy  in  drag  difference  is  due  to 
errors  in  the  estimated  incidence  at  a  given  lift.  As  noted  earlier,  the 
calculation  method  overestimates  lift  at  a  given  incidence  for  flows  with 
shock-induced  separation,  and  Wing  B  is  less  prone  to  shock-induced  separ¬ 
ation  chan  Wing  A.  Hence  there  appears  to  be  a  link  between  the  larger 
drag  difference  of  the  measurement  compared  with  that  of  calculation  and 
shock-induced  separation.  Therefore  it  is  reasonable  to  infer  that  shock- 
induced  separation  should  be  prevented  if  the  aim  is  to  minimise  drag.  It 
will  be  recalled  from  section  4  that  this  is  achieved  if  the  Mach  number  of 
flow  just  upstream  of  and  normal  to  the  shock  is  less  than  about  1.3. 

Finally,  Fig  21,  shows  the  result  of  an  analysis  of  the  lift-dependent 
drag  of  Che  two  wings  derived  from  the  experimental  results.  The  data  are 
presented  in  the  form  of  a  figure  of  merit  X  ,  the  ratio  of  the  lift- 
dependent  drag**  of  the  notional  wing  alone  to  that  of  linear  theory.*** 

The  calculations  by  linear  theory  are  made  for  a  flat  wing  of  cropped  delta 
planform  (Ref  18),  approximating  to  the  planform  of  the  present  wings,  as 
shown  in  Fig  21.  Consistent  with  the  definition  of  notional  wing-alone 


*  In  the  calculations,  drag  is  determined  by  integrating  the  streamwlse 
components  of  pressures  around  the  wing  and  combining  this  with  an 
estimated  skin-friction  drag.  The  calculated  drag  is  therefore  sen¬ 
sitive  to  errors  in  pressure,  particularly  near  the  leading  edge. 

**  The  zero-lift  drag  used  to  calculate  lift-dependent  drag  was  determined 
from  the  tests  on  the  third,  symmetrical  wing  of  the  family  of  wings. 

***  The  ratio  X  may  be  considered  to  be  a  figure  of  merit  since,  in  the 
linear  theory  approximation,  effects  which  reduce  the  thrust  forces  on 
the  wing  (eg  upper-surface  shock  waves)  arc  absent. 
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forces,  only  forces  on  the  part  of  the  vlng  aft  of  the  junction  of  the  wing 
leading-edge  and  the  body  are  considered  in  the  linear-theory  calculations. 
Also  shown  in  Fig  21  is  a  curve  derived  from  information  in  Ref  6 
illustrating  a  'practical  goal'  for  the  ratio  X  .  As  might  be  expected 
from  Fig  20,  Wing  B  is  superior  to  Wing  A  insofar  as  the  former  wing  yields 
lower  values  of  X  than  the  latter.  However,  for  values  of  C^u  greater 
than  about  0.14,  Wing  B  also  has  a  higher  lift-dependent  drag  than  that 
suggested  as  a  practical  goal.  For  values  of  C^w  greater  than  about  0.25 
this  can  be  explained  in  part  by  the  appearance  of  shock-induced  separation 
while  for  0.15  <  CLw  <  0.25  the  extra  wave  drag  due  to  upper-surface  shocks 
would  appear  to  be  the  sole  cause.  According  to  an  unpublished  theoretical 
study  by  Miss  C.J.  Betts  (RAE  Bedford)  of  wave  drag  of  conical  wings  of 
differing  camber  and  thickness,  wave  drag  due  to  upper-surface  shocks  is 
significant  only  for  values  of  %  greater  than  1.2.  Calculations,  based 
on  measured  pressure  distributions,  indicate  that  Mg  is  greater  than  1.2 
for  stations  aft  of  SOX  chord  on  Wing  B  for  Clw  “  0.2.  Thus  it  is 
reasonable  to  expect  further  reductions  in  lift-dependent  drag  for  lift 
coefficients  of  interest  if  %  is  reduced  to  values  below  1.2. 


6  CONCLUSIONS 

A  wind-tunnel  investigation  has  been  described  concerned  with  the  study  of 
flows  at  supersonic  speeds  over  two  highly-swept  wings  of  differing  camber 
but  with  the  same  thickness  distribution  featuring  rounded  leading-edges. 
The  camber  of  one  wing  (A)  was  derived  from  that  of  a  wing  designed  by 
linear  theory  while  the  second  (B)  was  designed  using  the  full-potential 
method  COREL  on  a  'net  wing'  basis  in  the  absence  of  a  representation  of  a 
body  in  the  method.  The  main  conclusions  drawn  from  the  study  are  as 
follows: 

(1)  The  approach  to  the  design  of  Wing  B  resulted  lr.  an  unexpected  re- 
expansion  on  the  upper  surface,  leading  to  shocks  which  were  strong 
enough  to  separate  the  boundary  layer  at  a  typical  supersonic 
manoeuvre  condition. 

(2)  Flow  separation  from  the  leading  edge  was  avoided  on  both  wings  at  a 
typical  manoeuvre  condition.  At  low  lift,  lower-surface  separation 
near  the  leading  edge  was  absent  on  Wing  B  but  was  Indicated  on  Wing  A 
by  a  collapse  in  the  lower-surface  suction  peak. 

(3)  For  Mach  numbers  of  the  local  f low-component  normal  to  and  Just 
upstream  of  the  upper-surface  shock  (Mg)  greater  than  1.3,  shock- 
induced  separation  was  observed.  Comparisons  between  measurement  and 
calculation  suggest  that  such  separations  on  the  upper  surface  Induce 
a  loss  of  lift  at  a  given  Incidence  on  the  part  of  the  wing  inboard  of 
the  shock,  causing  an  Increase  in  drag  at  a  given  lift. 

(4)  The  lift-dependent  drag  of  the  newly-designed  Wing  B  is  lower  than 
that  of  wing  A  and  further  reductions  in  drag  may  be  achieved  in 
future  wing  designs  by  ensuring  that  Mg  is  less  than  1.2. 

(5)  The  ARA/BAa  multiblock  system  for  solving  the  Euler  equations  has  been 
validated  for  both  wings  and  has  been-shovu  to  provide  a  sound  basis 
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for  designing  wing-body  combinations  for  efficient  manoeuvre  at  super¬ 
sonic  speeds  provided  account  is  taken  of  likely  boundary-layer 
effects. 
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Fig  5  General  layout  of  model  showing  pressure  plotting  locations 
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Fig  10  Flow  visualisations  and  upper  surface  isobars.  Wing  A, 


Fig  12  Flow  visualisations  and  upper  surface  isobars,  wing  B, 


Fig  14  Flow  visualisations  and  upper  surface  isobars.  Wing  B, 
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Fig  17  Comparison  of  spanwise  pressure  distributions 
Wing  B,  M  -  1.61,  Alpha  *  7.45 
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Fig  21  Lift-dependent  drag  against  lift  -  a  figure  of  merit  approach 
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